Quantum Tunneling and Relaxation in Mni2ac Studied by Magnetic Spectroscopy 
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The spin relaxation in the molecular magnet Mni2 ac is investigated by a novel type of high- 
frequency magnetic spectroscopy. By reversing the external magnetic field the intensities of the 
I ± 10) ^ I ± 9) transitions within the ground multiplet redistribute; below T « 2 K resonant 
quantum tunneling is directly observed as magnetic hole burning in the absorption spectra of a 
polycrystal. We quantitatively describe our observations by taking into account inhomogeneous 
line broadening, the time dependence of the level population during the relaxation process, and the 
resonance behavior of the relaxation rate due to quantum tunneling near the avoided level crossing. 

PACS numbers: 75.50.Xx, 78.20.Ls, 75.45.+j, 71.70.-d, 76.30.-v 
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In the last years molecular magnets comprising large 
clusters of coupled magnetic ions have attracted con- 
siderable interest as mesoscopic systems which exhibit 
new properties such as molecular magnetic bistability, 
macroscopic quantum tunneling of magnetization, quan- 
tum phase interference, etc. Among them clusters 
with a high-spin ground state and large anisotropy, such 
as Mni20i2(CH3COO)i6(H20)4]2CH3COOH-4H20 (ab- 
breviated as Mni2ac) are of a paramount interest. 

Mni2ac cluster can be considered as nanoparticle with 
an effective magnetic moment of {hb is the Bohr 

magneton), corresponding to a collective spin S" = 10 of 
the exchange coupled twelve Mn ions within the molecule. 
Due to the strong anisotropy, the magnetic moment 
freezes along one of the two easy directions at low tem- 
peratures. This is governed by a crystal field (CF) which 
splits the ground 5 = 10 multiplet and provides a quasi- 
doublet structure | ± m) of the lowest energy levels. The 
energy barrier /S.E « 60 K between the two lowest ly- 
ing states I ± 10) determines the thermally activated re- 
laxation of the magnetization in Mni2ac [^J]. Below a 
blocking temperature of about 3 K resonant quantum 
tunneling was discovered in the form of peaks in the re- 
laxation rate and steps in the hysteresis loops at regular 
intervals of magnetic field corresponding to a coincidence 
of the CF energy levels [|l],|j . Since the rate of the tun- 
neling for a ground state is small in Mni2 ac the tunnel- 
ing occurs via appropriate thermally excited CF states 
for temperatures above approximately 2 K resulting in 
a shortcut at the top of the barrier which accelerates 
the relaxation; at even lower temperatures pure quan- 
tum tunneling is dominant 



Various spectroscopic investigations have been applied 
to Mni2 ac, including electron-spin resonance (ESR) 
and inelastic neutron scattering [Q. In particular, a new 
type of high-frequency ESR method based on a quasi- 
optical spectroscopy Q made it possible to directly ob- 
serve the CF transitions in an equilibrium state and to 
study in details their characteristics, including the line- 



shapes ^ . In this Letter we report on the high-frequency 
magnetic spectroscopy of non- equilibrium phenomena in 
Mni2ac which allows to directly and independently de- 
termine the full set of magnetic characteristics. 

Mni2ac was pressed into a pellet with a thickness of 
1.39 mm and a diameter of 10 mm. The high-frequency 
ESR experiments were performed by means of a quasi- 
optical coherent-source technique . The transmission 
spectra Tr(j/) of the pellet were measured with a linearly 
polarized radiation (its magnetic field denoted as h), in 
the frequency range = 8 to 12 cm^^ at temperatures 
down to 1.8 K and in an external magnetic field H up 
to 8 T applied perpendicular to the radiation k-vector 
(Voigt geometry) . One frequency scan over the full range 
takes about 20 seconds. 

Examples of the equilibrium Tr(j/) spectra are shown in 
Fig. I^a. The absorption line at j/q = 10 cm~^ corresponds 
to the I ± 10) ^ I ± 9) CF transitions at H=0. It shifts 
to higher frequency VHiH) w vq + j\H,,\ (7 = giJ-s/h, 
Hz = H cos 9h, Oh is the angle between H and the easy 
C4 z-axis, g ~ 2) and changes its shape with increasing 
H, due to the inhomogeneous splitting of the degenerate 
I ± 10) and I ± 9) states in the polycrystalline Mni2 ac. 

To study metastable states of Mni2 ac we have investi- 
gated the time dependence of the Tr{i/) spectra. Apply- 
ing an external field H = +0.45 T shifts the absorption 
line to the position uh > 1^0 (Fig. |l|b); we then reverse 
the field (in about a minute), which leads to the line shift 
to the position ^ i/q ~ l\Hz\ < i^o- Starting from this 
point (t — 0) the intensity of the ^L-mode decreases with 
time while that of the i^/f-mode recovers correspondingly. 

The energy levels of Mni2 ac are determined by the ef- 
fective Hamiltonian (|-|^ 

H,s = DSl + BSt + C{Sl + St)/2 - g^sS ■ H (1) 

where S is the effective spin of the ground S* = 10 mul- 
tiplet, and D, B and C are the CF parameters. At small 
fields the low energy levels are determined mainly by Hz '■ 
Em — Dm^ + Bm'^ — gfisHzm, m = 0, ±1 . . . ± 10 (inset 
of Fig. lb). At low enough temperatures for Hz > only 
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the I + 10) quantum state is populated and we observe 
only one magnetic dipole transition | + 10)^| + 9). After 
reversing the field {Hz < 0) the state | + 10) becomes 
metastable and its population relaxes to the new ground 
state I — 10); the intensity of the transition | + 10)^| +9) 
(at < i^o) decreases while that for the transition 
I — 10)— >| —9) (at vh > t'o) grows (Fig. lb). For our poly- 
crystalline non-oriented sample a random orientation of 
the easy-axis results in the population of both | -I- 10) and 
I — 10) CF states and in an inhomogeneous broadening 
of the corresponding transitions by the applied field. 

New features in the relaxation phenomena appear at 
lower temperatures, T < 2 K, as shown in Fig. ^. The 
relaxation occurs "inhomogeneously" ; i.e., it has a faster 
rate near a particular frequency i^i^^ = 9.57 cm~^ where 
a peak in the Tr(t^) spectrum appears accompanied by 
a minimum developing at = 10.47 cm~^. This new 
phenomenon can be called magnetic hole burning in anal- 
ogy to spectral hole burning in dielectrics [|ll|. We as- 
cribe this fast relaxation near the v^^"^ to a resonant quan- 
tum tunneling in those crystallites whose orientation rel- 
ative to H satisfy the level crossing condition hu'^\m) — 
hvo ± gnBHf°^ = hvo ± [-D - B{w? + {m - if)]. 

For a quantitative description of the observed relax- 
ation the effective magnetic permeability is calculated 
using (a) the non-equilibrium populations of the low CF 
states I ± 10) and their evolution (relaxation) with time 
and (b) the inhomogeneous broadening of the CF transi- 
tions in the magnetic field due to the random orientation 
of the crystallites. Firstly, we analyze a behavior for a 
single crystallite. Considering only two | ± 10) ^| ± 9) CF 
transitions and introducing a non-equilibrium normalized 
populations Pm{t) we generalize the ordinary equilibrium 
permeability and obtain its transverse components: 

Hxx.,yy{v,t) = fj.±{iy,t) = 1 + Afi'^{i^,t) + /\^^{v,t) (2a) 

IJ.xy,yx{v,t) = ±Ha{v,t) = ±i [A/i+ (z^, t) - A/i" (l^, t)] (2b) 



where 



A^±(iy,i) = A^o 



(3) 



are the contributions of the | ± 10) ^| ± 9) CF transi- 
tions to the permeability with its static value at T = 0, 
H = being Apo; a Lorentzian lineshape R^{v) is as- 
sumed for simplicity; hv± = E±q — E±iq — hvQzLgpBHz- 
The equilibrium state is given by Prn.{t oo) = p"^ = 
exp{-PEm}/Z with Z = J2eM~PEra}, 13 = l/fcsT, 
and Y^Pmit) = 1. 

In the thermally assisted regime, the relaxation of the 
populations Pm{t) takes place via excited states and is 
determined by the complex master equation for the pop- 
ulations of all CF states |12-1^]. It can be reduced to 
the simple relaxation equation Ap — ~{Ap — Ap°°)/T 



difference of the equilibrium population after the field is 
reversed. The solution Ap{t) — Ap°°[l - 2exp{— t/r}] 
for the initial condition Ap{Q) — Ap° — —Ap°° leads to 
the population differences 

Ap±{t) = p±io - p±9 = Apt + {Ap\ - Apt) e-*/^ (4) 

which enter Eq. (||); here Apt ~ [exp{/3 i?±io} — cxp{- 
/3£;±9}]/Z and Ap\ = [1 T A/?°°(l + exp{-/3/ii/zp})] (1 - 
exp{— /3/ii/±})/[2 -I- exp{— /3/ii/+} -|- exp{— /3/iz^_}] , where 
all quantities refer to the reversed field. The deriva- 
tion assumes thermal equilibrium between | ± 10) and 
± 9) states due to fast relaxation via phonons, im- 
plying p±g = p±io exp{— For low temperatures 
{f3hh'± >> 1), p±Q is negligible and Eq. (^) simplifies to 
Ap±{t) « p±io(i) « ptw + {Ptw - P±w) exp{-t/T}. 

All previous consideration refers to a single crystal. 
For polycrystalline samples we have to take into ac- 
count the inhomogeneous broadening of the lines due 
to different splittings of the doublets. By averag- 
ing the permeability [Eq. (||)] of the crystallites with 
respect to their orientations we obtain the effective 
permeability p'^^{v,t) with the non-zero components 

Pfx^yy{^,t) = {p±{iy,m + COS^0h))/2, p%%^yA^.t) = 

±i(^a(i/, i)cos0//), and p°^{v,t) = {p_L{v,t)s\v? Oh) for 
the H oriented along the z-axis. Using pp^ [v, t) and the 
Fresnel formulas for the transmission coefficient jl^ we 
fit the zero-field spectrum and obtain vq — 10.02 cm~^, 
F — 0.18 cm^^, and Apo = 0.0107, in accordance with 
our previous results These parameters also describe 
the spectra in the magnetic field |lj] in all details, in- 
cluding the dramatic changes of the lineshape (Fig. |^). 

The relaxation phenomena observed in our Tr(i^) spec- 
tra at various temperatures can be self-consistently de- 
scribed using the theory of the magnetization relaxation 
in Mni2 ac |l2|-p^ based on phonon-assisted spin tunnel- 
ing induced by forth order CF and transverse magnetic 
field. To determine the relaxation rate t{T, HjCOsOh) 
we have numerically diagonalized the master equation for 
the density matrix p = Wp in the space of the ground 
S — 10 multiplet, following the procedure suggested 
by Ref. ||l^. In calculations we used the CF parame- 
ters of Eq.(l) D =-0.389 cm^i, B =-7.65x10"'' cm-\ 
C| = 5.0 X 10^^ cm^^ from (|-||] and parameters of 
the spin-phonon interaction determined via the main CF 
term DS^ 
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In Fig. ^ the angular dependence of the relaxation time 
is plotted for two different temperatures and fields used 
in our experiments (Figs. ^ and ^). The relaxation goes 
fastest for cos 9h — because the energy levels on the 
both sides of the barrier coincide; a situation, however, 
not observed in Voigt geometry. For fields exceeding the 
value of the first level crossing (~ 0.44 T), additional 
resonance minima arise with a fine structure related to 
the slightly different resonance conditions for various ex- 
where Ap°° — pfio ^ P^w is the citing levels due to the fourth order CF contribution B. 
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As an example, Fig. |i|b shows the relaxation observed in 
the Tt{v) spectra at iJ = 0.45 T and T = 2.6 K which 
are well described by our model. We note that a weak 
resonance near cosOh = 0.97 due to avoided level cross- 
ing of the CF states to = 4, —3 (or to = 3, —4) does 
not exhibit in the spectra due to the averaging in the 
polycrystalline sample. The relaxation rate strongly in- 
creases with magnetic field because the effective barrier 
decreases with both (longitudinal and transverse) com- 
ponents 1^; there are also significant contributions from 
the tunneling processes which are enhanced because the 
tunneling splitting increases with the transverse field. 

These resonance effects are best seen at low temper- 
atures, when the thermal relaxation over a barrier is 
suppressed. Fig. ^ shows a set of resonances in the re- 
laxation rate at T =1.96 K which correspond to pairs 
of states with an avoided level crossing. Using this 
t{cos9h), the time evolution of the Tr(j/) spectra is 
calculated; as seen from Fig. || the agreement is excel- 
lent. The main contribution to the resonance relax- 
ation comes from the thermally assisted tunneling via 
(5,-4), (6,-5) and in less degree (7,-6) states (or cor- 
responding states with the opposite signs); they deter- 
mine the frequency of the burned spectral hole. For 
these avoided level crossings the resonance lineshape is 
determined by a truncated Lorentzian [ p^JT5| with the 
effective linewidth determined by the tunneling splitting 
(0.23, 0.015, and 5.0 x 10~^ cm~-^, respectively) which 
is much larger than the corresponding phonon linewidth 
(ranging from 5 x 10^^ to 3 x 10^^ cm^^). The linewidth 
(0.12 — 0.18 cm^^) of the zcro-ficld-resonance absorption 
vq exceeds the corresponding phonon linewidths consid- 
erably. It is mainly determined by inhomogeneous broad- 
ening and thus serves as a measure of the disorder effects 
1^ , which could originate from CF and g- factor distribu- 
tions (D-strain and g-strain) and random magnetic dipo- 
lar fields as was shown by ESR 

The problem of shape of the tunneling resonance 
and the effect of inhomogeneous broadening on it are 
presently subject of considerable interest and discussions 
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I 19|-|2^, many issues, however, remain open. One 
of them could be related to the obvious discrepancy be- 
tween the observed inhomogeneous character of the zero- 
field modes broadening and Lorentzian shape of the 
tunneling resonances observed at temperatures T > 2K 
|^ , pT| and described theoretically in the thermally as- 
sisted regime [Q. One may expect that a large inho- 
mogeneous broadening of the levels (i.e., larger than the 
tunneling splitting) reduces the tunneling amplitudes and 
smears of the resonances, resulting in suppression of the 
relaxation peaks and in change of their Lorentzian shape 
I psf . Suggesting that the main contribution to the line 
broadening comes from the CF dispersion and assuming 
its Gaussian distribution [ p2[ , we can evaluate the dis- 
persion 6D — 0.004 cm^^of the main quadratic term D 
in the Hamltonian d|) [E5| . This yields a smearing of the 



tunneling level differences (5(i?m — £'-m+i) = 5D(2m—l), 
for instance, of 0.035 — 0.05 cm~^ for (5,-4), (6,-5), and 
(7,-6) states. From a comparison with the calculated 
tunneling splittings we can conclude that the inhomo- 
geneous level broadening does not effect the main reso- 
nance (5,-4), smears the (6,-5), but suppresses the (7,-6) 
and higher order resonances. A numerical simulation in- 
cluding an additional averaging with respect to the CF 
distribution confirms this conclusion. Thus, if the tun- 
neling splitting is sufficiently large for states contributing 
to the thermally assisted tunneling, the lineshape of the 
corresponding resonance remains Lorentzian. 

In summary, we utilized a novel optical method of 
high frequency magnetic spectrosopy to study the re- 
laxation phenomena in Mni2ac clusters, including the 
resonance quantum tunneling along with the CF tran- 
sitions and their lineshapes and linewidths. At elevated 
temperatures only frequency-homogenous relaxation due 
to mainly thermal activation is observed, whereas for 
T < 2 K the thermal relaxation is significantly sup- 
pressed and resonant quantum tunneling dominates; in 
this range the effect of magnetic hole burning is observed. 
We have developed a microscopic model which allows to 
quantitatively describe all observed effects, taking into 
account the orientational distribution of the relaxation 
time and to determine the most effective relaxation chan- 
nels and to estimate the effect of CF inhomogeneity on 
the resonant tunneling. 
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FIG. 1. Transmission spectra of poly crystalline 1.39 mm 
thick pellet of Mni2ac measured in an external magnetic field 
H ± h _L k. (a) Shift of the absorption line with increas- 
ing magnetic field. The solid lines represent Lorentzian fits 
with the distribution of the crystallites orientation taken into 
account, (b) Time evolution of the absorption line after the 
magnetic field is reversed from -1-0.45 T to —0.45 T. The lines 
represents the transmission for the calculated angular distri- 
bution of the relaxation time t{cos6h) exhibit in Fig. 3. The 
arrows indicate the positions of line at H=0, ±0.45 T. Inset: 
the energy levels of a Mni2 ac single crystal for a magnetic 
field H II C4; the arrows show the | ± 10) — > | ± 9) transitions 
before and after the field inversion. 

FIG. 2. Time evolution of the transmission spectra at 
T — 1.96 K after the magnetic field (H _L h ± k) is reversed 
from -1-0.9 T to —0.9 T. The narrow absorption dip within 
the broad feature is due to resonant quantum tunneling be- 
tween pairs of coinciding energy levels. The lines correspond 
to the transmission spectra for calculated distribution of the 
relaxation time t{cos6h) shown in Fig.3. Also shown is the 
zero-field absorption line. 

FIG. 3. Calculated distribution of the phonon-assisted 
spin-tunneling relaxation time t{cos9h) as a function of the 
angle between the magnetic field H and the C4 axis of the 
crystallites at T = 2.6 K and T = 1.96 K for certain values 
of the applied magnetic field. The numbers at the resonance 
minima indicate pairs of tunneling CF states. 
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